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1 Introduction 

In the past, the structure and its control system have been designed independently. Structural design and 
optimization, and control system design and optimization, have each been areas of separate research, each 
progressing vigorously along its own path. However, spurred on by recent proposals of new, large, highly 
constrained space structures, the question has arisen as to whether an integrated structural/control design 
procedure might not be more appropriate. The first papers actively investigating the question of simultaneous 
structure and control design began appearing in the literature around 1983 [1,2,3]. Since then, there has been 
a growing interest in this subject from other authors, although the field itself is still in relative infancy. Using 
a conventional design approach for a controlled structure, one would first optimize the structure alone, then 
design a control system for this baseline structure. This process may then be iterated until both the structure 
and control system meet necessary constraints and objectives. 

Some authors ([4, 5, 6, 7], for example) take a “classical” approach to the simultaneous structure/control 
optimization by attempting to simultaneously minimize the weighted sum of the total mass and a quadratic 
form, subject to all of the structural and control constraints. In this paper, the optimization will be based 
on the dynamic response of a structure to an external unknown stochastic disturbance environment [8]. Such 
a “response to excitation approach” is common to both the structural and control design phases, and hence 
represents a more natural control/structure optimization strategy than relying on artificial and vague control 
penalties. The design objective is to find the structure and controller of minimum mass such that all the 
prescribed constraints are satisfied. 

Two alternative solution algorithms will be presented which have been applied to this problem. Each 
algorithm handles the optimization strategy and the imposition of the nonlinear constraints in a different 
manner. Two controller methodologies, and their effect on the solution algorithm, will be considered. These 
are full state feedback and direct output feedback, although the problem formulation is not restricted solely to 
these forms of controller. In fact, although full state feedback is a popular choice among researchers in this field 
(for resons that will become apparent), its practical application is severely limited. The controller/structure 
interaction is inserted by the imposition of appropriate closed-loop constraints, such as closed-loop output 
response and control effort constraints. Numerical results will be obtained for a representative flexible structure 
model to illustrate the effectiveness of the solution algorithms. 

2 General Problem Formulation 

The integrated control/structure design optimization problem can be stated as follows: find the vector of 
structural and controller parameters that minimizes the mass of the structure subject to a set of prescribed 
stochastic disturbances, with limitations on the available control energy and on a set of allowable output 
responses. This can be written in the form of a nonlinear mathematical programming problem as 
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is an iV-vector of design variables, 

is an m-vector of structural constraints, 

is an n- vector of state variables, 

is an n u -vector of control forces, 

is an n w -vector of stochastic disturbances, 

is the (n x n) matrix containing the system dynamics, 

is the (n x n u ) matrix containing information on the 

locations and orientations of the actuators, 

is the (n x n w ) matrix containing information on the points 

of application and orientation of the disturbances, 

is the i th control effort constraint cost function, 

is the n u .-order partition of u representing the control 

forces involved in </ cc ., 

is the maximum allowable value of the i th expected control 
effort function E[tif 

is an (n U| x n u .) control force weighting matrix, 

is the i ih output response constraint cost function, 

is the maximum allowable value of the i th expected output 

response function E [y£ .], 

is an (n</ . x n<*. ) output response weighting matrix, 

is the i th design output -vector, 

is an (n<*. x n) matrix giving the relationship between 

the state variables and y d . , 

is an JV-vector of minimum design variable values, and 
is an N-vector of maximum design variable values. 


The side constraints are the strict bounds p t and p u on the design variables, and are vector inequalities 
that are imposed element by element. These design variable bounds are not included explicitly as constraints 
in the problem formulation. Note that the structural weight and the structural constraints g will in general 
not be functions of the controller design variables unless the controller mass is included in the design. Note 
also that g CCi is a weighted mean square control effort, and g 0Ci is a weighted mean square output response. 
Multiple output response constraints are allowed, although only one of these will in general be active at the 
optimum design. However, all of the control effort constraints will generally be active at the optimum design. 

In this work w is a zero mean Gaussian white noise disturbance with covariance X w . The structure will 
respond to this disturbance with some transient behaviour, in addition to a steady-state response. It seems 
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reasonable to optimize the structure for the steady-state response to the disturbance rather than the transient 
response because the transient behaviour will normally be of secondary importance to the response objectives 
(such as long term pointing accuracy). Also, for steady state optimization, the differential equation constraint 
(state equation) can be replaced with a steady state covariance equation, so that the control effort and output 
response constraints may be recast in terms of this covariance. Therefore the two-point boundary value problem 
is eliminated and the numerical solution of the problem is significantly simplified. 

2.1 Full State Feedback Control 

The simplest form of feedback control is to feedback the entire state vector, with u — —Kx, where K is the 
(n u x n) state feedback gain matrix. The controller design variables for this case will be the n u n elements 
of K. Substituting this control into the state equation, and assuming that the disturbance w is zero mean 
Gaussian white noise, the state covariance matrix X for this case can be found from the Lyapunov equation 

F el X 4- XF% + G w X w Gl = 0 (2) 

where F c i = (F — GK) is the stable closed-loop dynamical matrix for the full state feedback case, X = E[xx T ] 
is the (n x n) symmetric state covariance matrix, and X w = E[tuio T ] is the ( n w x n w ) symmetric covariance 
matrix for the stochastic disturbances. 

Expressions for the controller constraints in terms of this covariance matrix can then be obtained as 
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where Ki is the ( n Ui x n) partition of K corresponding to It is assumed that the ti, are independent, and 
that u and K are ordered as 


« T = [«T ul ul' ], K T = [Kf Kl ••• /C, ] (5) 

Note that Ya=i n »i = n «> an <i that the columns of G can be interchanged to force condition (5) to be satisfied. 
Using full state feedback, the first-order necessary (Kuhn-Tucker) conditions for optimality can be analytically 
solved to give [8] 


K = R~ l G T A^, where F T A X -f A X F — A X GR l G T A x + W = 0 (6) 

and where R and W are respectively the (n u x n u ) and (n x n) matrices defined as 

The variables A u - and A y , come from the Kuhn-Tucker conditions, and are the Lagrange multipliers associated 
with the i th control effort and output response constraints respectively. 

Equations (6) define the solution to the optimal control problem 

min J c = / [* T Wx + u T Ru ] dt (8) 

Jo 

where K is the optimal steady-state gain matrix, and A x is the steady-state solution to the associated Riccati 
equation. Although this LQR property only holds true at the optimum point, it is computationally convenient 
to assume that at every point in the design cycle, the control design variables will be found as the solution to the 
optimal control problem (8). Therefore, the numerical optimization problem can be reduced to optimization 
over just the structural design variables, along with an optimal control problem solution which will be a function 
of the Lagrange multiplier vectors A u and A y . The immediate benefit of this is a reduced dimensionality 
nonlinear programming problem. In addition, since the regulator solutions always give a stable closed-loop 
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system, no explicit check must be performed on the system stability during the solution procedure. The LQR 
assumption in this proble formulation is similar to the approach taken in [9,10], and others, where R and W 
are fixed, and not chosen to satisfy the constraints. 

2.2 Direct Output Feedback Control 

For most real systems, the state vector will be very large, and the use of full state feedback would result in a 
controller of unacceptably high dimension, assuming additionally that the entire state is available. However, 
usually only a small subset of the system states will be available to the designer, in the form of the output 
measurement vector. These can include actual system states along with linear combinations of the system 
states, in the form y=Hx , where y is the n y -ve ctor of outputs and H is the (n y x n) output matrix giving the 
relationship between the outputs and the system states. If these output states are to be used in the feedback 
loop, the resulting control is termed direct output feedback, with the control forces defined to be u = -Ky, 
where I< is the (n u x n y ) output feedback gain matrix. In this case, the controller design variables will be the 
n u n y elements of K . 

Substituting this control into the state equation, the state covariance matrix X for this case can be found 
as the solution to the same Lyapunov equation (2), except that now the closed-loop dynamics are given by 
Fci = (F — GKH). Note that since K does not satisfy any special conditions (such as the LQR conditions), 
the closed- loop system F -/ is not guaranteed to be stable for any K . If F e i is unstable at any stage in the 
solution procedure, the covariance matrix X cannot be found from equation (2). Therefore, for this case, hard 
constraints on the closed-loop system eigenvalues must be imposed at every step in the design proceedure. 
“Hard” in this sense means that special precautions must be taken in the solution procedure such that these 
constraints can never be violated. 

Expressions for the controller constraints for direct output feedback in terms of the covariance matrix can 
then be found to be 
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3 Solution Algorithms 

In the general problem formulation presented in the previous section, the constraint functions are generally 
highly nonlinear implicit functions of the design variables. Solution of this problem could be attempted by the 
direct application of nonlinear programming techniques; that is, using the exact functional expressions for the 
constraints. However, this approach quickly becomes computationally very expensive as the dimensionality 
increases since the full objective and constraint functions must be evaluated at every step, and their respective 
gradients at most, if not all, steps thoughout the design procedure. Such evaluations tend to be computationally 
very expensive. 

Approximation techniques, where the implicit nonlinear problem is replaced by a sequence of explicit 
approximate (although not necessarily linear) problems, have been shown to yield efficient and powerful al- 
gorithms for structural design optimization (see, for example, [11,12]). In this paper, two solution techniques 
based on approximation techniques will be tested on the integrated control/structure design optimization prob- 
lem. The methods will be compared with respect to the ease of use, generality of application, and numerical 
robustness to changes in move-limits and other solution parameters. 

3.1 Sequential Nonlinear Approximations 

In this method, the fully constrained nonlinear optimization problem is solved by the iterative construction 
and numerical solution of a sequence of explicit approximate problems. The approximate problems are first- 
order Taylor’s series expansions (with respect to either the inverse design variables ([13], for example), or with 
respect to hybrid design variables [14]) of the objective and constraint functions. Depending on the intermediate 
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variables chosen, the approximate functions may still be nonlinear functions of the design variables. Therefore, 
the numerical solution is accomplished using a mathematical programming code, specifically the modified 
method of feasible directions as implemented in ADS [15]. 

The solution process begins with some initial structure, which is analyzed using the finite element technique. 
At this point, the gradients of the active constraint set are evaluated, and the approximate problem is formed, 
with respect to the current design. Expressions for the gradients of all constraints considered can be evaluated 
analytically. The approximate problem is solved with ADS using an active constraint set strategy to reduce 
the dimensionality of the approximate problem by deleting the inactive constraints. Move-limits on the design 
variables are imposed during the solution to ensure that the design remains within the region for which the 
approximation functions are of acceptable quality. The choice of move-limits and how they change can have a 
significant effect on convergence, and will often be determined from numerical experience with the particular 
problem at hand. 

After the solution of the approximate problem, the structure and its control system are deemed optimal 
if a convergence test on either the absolute or relative objective function change over a specified number of 
successive global iterations is satisfied. Otherwise, the objective and active constraint gradients are evaluated 
for the new design, a new approximate problem formed, and the process above is repeated in an iterative 
manner. The solution procedure ends when the design variables converge, or when the number of iterations 
exceeds some preset maximum. 

Scaling the structure and controller to the closest constraint surface may be possible in some cases, because 
of the special assumed form of the controller. Scaling to structural constraints has been performed in other 
work (see [16]) and will not be covered here. If full state feedback is used, it is possible and practical to scale 
the structure to the closest control effort constraint and closest output response constraint simultaneously. The 
variables with which the structure is scaled are the structural design variables (elemental areas or thicknesses), 
and the Lagrange multipliers associated with the two controller constraints A u and A y (where for clarity, and 
without loss of generality, the subscripts on the A’s that refer to the particular control effort or output response 
constraints under consideration have been dropped). 

Note that changing the values of A u and A y cannot independently change the values of u m , = tr (K T RKX) 
and y m3 = tr(HjW HdX), because in the LQR problem, only the ratio of A u to A y is important. One can 
choose the ratio (A u /A y ) to satisfy one of the control constraints — say u mj> . Then y ms will not in general be 
satisfied. Suppose y ms is too large (i.e. y m8 > a 2 ) at the particular point where u ms is satisfied. Then the only 
way one can satisfy the y ms constraint is to increase the sizes of at least some of the structural members. This 
seems reasonable because if the control constraints could be satisfied by simply choosing appropriate controller 
parameters, then there would be no interaction between structural optimization and controller optimization. 
Intuitively, it can be seen that this is not the case. Note that each member of the structure will be scaled 
by the same amount to fulfill our goals. Obviously, this method is not absolutely mandated, and some other 
approach could be used where the design variables are not scaled equally. However, this would then be resizing 
rather than scaling , a process normally left to the nonlinear programming algorithm. 

The final scaling aim is to set u m8 = (3 2 and y ms = a 2 . To perform the scaling, it is assumed that, at 
iteration i, the values (t* m *)t and (t/m<)* will change, as a result of changes to (A u ),- and the (pj)i, according 
to the equations 
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If initial (educated) guesses for these constants can be made, they can be updated in an adaptive manner 
during the scaling proceedure. 

Move-limits are imposed on the design variables during each approximate problem solution. This is done 
in an attempt to restrain the design variables to a region in which the explicit function approximations remain 
reasonably accurate. However, deciding how to impose these move-limits is a non-trivial task. The local curva- 
ture of the design space (i.e. how nonlinear are the actual constraint surfaces in the region about the expansion 
point of the approximations) will determine the move- limits, with more strict move- limits applied in regions 
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of high curvature, and less strict move-limits imposed in regions of low curvature. Since second-derivative 
information is required to estimate curvatures, and since such evaluations are very expensive computationally, 
imposing move-limits is usually reduced to an art based on past experience. Quasi-Newton methods obtain the 
second derivatives using only first derivative information, however, these methods typically take N iterations 
to fill the Hessian, and can be very costly if N is large. 

For the purpose of this work, a move-limits factor 7 is imposed in an exponential form. If the current design 
variable and approximation expansion vector is p, then the upper and lower bounds on the design variables 
for the current approximate problem are defined as 


Pu = 7 P» 



(13) 


where 7 > 1. The limits specified in equation (13) must be imposed element by element. Note that since 
the design variables in this example will be structural design variables only, they are restricted to be positive. 
Obviously, equation (13) must be modified if the design variables can be negative. The exponential form of 
the move-limit factor is defined by the particular choice of 7 m i n and y maT (typically 1.2 and 1000 respectively 
in this work). 


3.2 Continuation and Sequential Linear Programming 

The complex nature of the constraint functions in the nonlinear optimization problem, especially the con- 
troller constraints, leads to various convergence problems in the context of a classical gradient based nonlinear 
programming code such as ADS. As the problem dimensionality increases, convergence will usualy become 
increasingly difficult to accomplish, as step sizes reduce to satisfy the local linearity assumptions inherent in 
gradient based solution techniques. Another method for the solution of mathematical programming problems 
that has recently become popular is the use of continuation methods to impose nonlinear constraints cou- 
pled with sequential linear programming (SLP) [17,18]. The continuation procedure is a conceptually simple 
method of applying restrictive constraints gradually from less restrictive ones, which replaces the most de- 
manding constraint functions of the form g(p) < 0, by a set of neighbouring constraint functions Qi , defined 
by 


Gi{p * . 7 i) = g(Pi) - (1 - 7i)g(Po) <0 for * = 0, . . . , M (14) 

where p 0 is the arbitrarily chosen' initial design point, and 7 * is a continuation parameter satisfying 

0 = 7o<7i < — <7„ = 1 (15) 

Note that for 70 = 0 , when p = p 0 , the new constraint function Q 0 is identically satisfied. If convergence 
is acheived for 7 = 1 , then the original constraints will be recovered in M steps. The step size A7 = 7,- — 7,-_i 
(and hence Af), can be chosen small enough so that assumptions on local linearity can be almost arbitrarily 
satisfied. 

Linear Programming (LP) methods are a powerful approach to handling a large number of locally linear 
constraints, and due to the wide availability of very efficient LP codes, are an attractive alternative to nonlinear 
programming methods. The neighbouring problems generated by the continuation procedure can be written 
in the form 


Minimize J(pJ, subject to Gi{Pi,*fi) < 0 (16) 

To transform these equations into a linear programming problem, the equations are linearized about the current 
point p t , and move-limits on the maximum parameter changes allowable locally are imposed. Expanding the 
objective and constraint equations in (16) to first order in a Taylor’s series expansion about p { gives the locally 
linearized problems in linear programming form as 
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Minimize, with respect to A p t -, A J{ = 


A p t , subject to 
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—c < A Pi < t 


(17) 


for i = 0, . . . , (Af — 1). All elements of c are assumed positive, and the vector inequality is imposed element 
by element. 

The algorithm begins with an initial structure, which is analyzed using the finite element technique. The 
initial problem (70 = 0) is solved, and the continuation parameter 7 is incremented from 7 0 = 0 by A7 to 
71. Note that if initially K = 0, so that the control effort allowed for the initial local problem is zero, this 
initial problem becomes a pure structural optimization subject to the dynamic output repsonse constraints. 
The increment A7 is set by an a priori choice of Af, the number of continuation steps, although A7 need not 
be constant throughout the solution procedure. Successful implementation of the continuation method has 
been reported when A7 was chosen as initially quite small and increased to a larger value during the solution 
[18]. The choice of A7 is closely coupled with the choice of the nominal design variable move limit vector £0* 
There is in fact a tradeoff between the satisfaction of the local linearity assumption through c, and the ability 
to converge to the neighbouring problem through A7. Usually, for each particular problem, some numerical 
trial and error will be required to find those values of Co and A7 that yield an efficient solution technique. 

The gradients of the objective and constraint functions are calculated at the current design, and then the 
associated linear programming problem (17) is solved by a linear programming code. In this work, the linear 
programming routine E04MBF from NAGLIB (National Algorithms Group LIBrary) was used, although other 
routines inserted at this point should provide the same solution. Since the local linearity assumption will 
never be exactly satisfied, the actual constraint values at the new point, specified by the solution to the linear 
programming problem, will be different than that predicted. Therefore, the constraints Qi may not be satisfied 
following the linear programming solution step. 

Since a converged subproblem solution is required before increasing the continuation parameter, this local 
problem is iterated locally until convergence is obtained. At each local iteration, new gradients are calculated, 
and the move limits on the design variables are reduced so that e = ce 0, where a value of c = 0.75 was used 
in this work. If convergence to the local problem does not occur within 15 iterations (where move limits are 
about 1% of their nominal values), the move limits are increased to their nominal values £o> and the local 
iterations are repeated. Numerical experience with this algorithm has shown that this procedure is flexible 
enough numerically so that converged subproblem solutions can be obtained in a reasonable number of local 
iterations, as long as the neighbouring problems are “close enough”. Practically, this means that either the 
constraint values of the initial system should be “close” to their final desired values, or that AT should be 
large. Once the local problem has been solved, the continuation parameter 7 is incremented, and the new local 
problem solved as before. At the M th continuation step, 7 = 1 and the original problem is recovered, so that 
the solution to the M ih local problem is the solution to the original problem. 

If closed-loop stability constraints are violated at any stage in the solution procedure, these must be 
imposed immediately. To achieve this, it is possible to employ a method that never requires the calculation of 
the closed-loop eigenvalue derivatives, saving considerable computational expense. Of course, if in addition to 
overall stability there are constraints on closed-loop damping ratios or bandwidth, then the evaluation of the 
closed-loop eigenvalue derivative may be necessary at some point. The method used in this paper is to simply 
bisect A p i and perform another analysis until a stable system configuration is obtained. 


4 Gradient Analysis 

For the numerical optimization procedure to be practical, especially as the dimensionality increases to realistic 
structures, it is essential that it be possible to evaluate the first-order sensitivities of the complex constraint 
functions in an efficient manner. The objective function (the weight) is a linear function of the finite element 
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thicknesses and/or cross sectional areas (for truss type finite elements), so that its gradient is easy to calculate 
at any point in the design space. 

4.1 Gradients for the case of Full State Feedback Control 

The gradients of the controller constraints with respect to a structural design variable pj are given by 
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Note that gradients with respect to controller design variables need not be evaluated since these design variables 
were effectively removed from the optimization problem by the LQR constraint. 

4.2 Gradients for the case of Direct Output Feedback Control 

The gradients of the controller constraints with respect to a structural design variable pj are given by 
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where y>i, Zi, and Nj are evaluated using the following set of equations: 
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The gradients of the controller constraints with respect to the elements of the gain matrix K can be written 
in matrix expression form as 
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Note that the Ri in equation (35) is of order (n Ui x n u .), and is in the i th diagonal block of lli. 
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5 Example: The DRAPER I Teterahedral Truss Structure 

The DRAPER I structure [19] is a tetrahedral truss attached to the ground by three right-angled bipods, as 
shown in Figure 1. Although attached to the ground, this model will act as a typical flexible structure pointing 
subsystem (e.g. antenna, radar, optical) attached to a rigid core. Any motion would then be with respect to 
this rigid core, and transmit forces to it. Consequently, this model has no rigid body degrees of freedom. The 
finite element model has 12 truss elements, since the joints are pinned and transmit no moments. There are 
four nodes that are free to move in all directions, so the model contains 12 degrees of freedom. The structural 
design variables are the cross-sectional areas of each of the 12 truss elements. Since there are 12 degrees of 
freedom in the model for this structure, the state-space model will be 24 th order. There will be 6 inputs 
corresponding to the 6 legs of the structure, and a varying number of outputs, depending on the problem at 
hand. 

For the purposes of this work, material parameters of p = 0.1 lb/in and E — Young’s Modulus = 20 kpsi 
were used. The dimensional values E and p were chosen to give initial numerical values of structural frequencies 
for the dimensional model roughly comparable to those of the non-dimensional model. The model contains 
no nonstructural mass. Elements 7 through 12, the three right-angled bipods, take on the duties of force 
actuators (and possibly colocated velocity and/or displacement sensors). Only one output response constraint 
is defined (n a = 1), with the design output vector y d representing the line-of-sight error of the top vertex [(z, y) 
displacements of vertex 1], The disturbances, labelled u>i and it >2 in Figure 2, are assumed to be independent, 
zero mean, Gaussian disturbances with intensity 1.0. 

The damping added to the state space system will depend on the state space realization used. For cases 
Where a realization based on physical variables is used, the damping matrix C is formed to be C = 0.1 M + 
0.001 K s , For cases where a realization based on modal variables was used, the damping ratio of each mode 
was specified to be 0.1% of the modal frequencies during the formation of the state matrices. The weighting 
matrices R and W are set to the identity matrices, so that equal weighting is given to all components of u and 
y d . The minimum cross-sectional areas for all elements was specified as 0.1 in 2 . For this problem, no static 
structural constraints were specified in this model of the DRAPER I structure, the intent being to investigate 
the effect of the closed-loop controller constraints on the structural design optimization. 


5.1 Full State Feedback Control 

In this section, the sequential nonlinear approximations solution algorithm, with the addition of the scaling 
procedure outlined in Section 3.1, is applied to the full state feedback control of the DRAPER I structure. 
The effect of the scaling procedure used here can be determined by applying the sequential approximations 
algorithm in the form of a direct output feedback problem with H - I with no scaling assumed. The 
continuation solution algorithm is also best handled in the form of a direct output feedback problem since 
no special scaling is assumed. Both solution algorithms applied to the case of full state feedback are discussed 
in Section 5.2, where direct output feedback control is considered. For brevity, only limited results for both 
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controller methodologies are presented in this paper, but a full discussion of this example can be found in 
reference [20]. 

Runs were made optimizing the DRAPER I structure using an inverse design variable approximation for 
all constraint functions. The initial structure was defined with all structural design variables set at 10 in 2 , and 
with the Lagrange multipliers A u and X y set at 1.0. This set of initial conditions will be termed the symmetric 
set of initial conditions, for they specify a structure with a number of vibrational modes of the same frequency 
(repeated eigenvalues). A range of allowable expected output response (a 2 ) of 1 x 10” 5 in 2 to 1 x 10“ 4 in 2 in 
steps of 1 x 10 -4 in 2 , and allowable expected control effort (/? 2 ) of 50 lb 2 to 80 lb 2 in steps of 10 lb 2 were used. 

Table 1 summarizes the resulting minimum weight in pounds found for the case where a state- space real- 
ization based on the modal displacements and velocities was used. Intuitively, two trends would be expected 
in the data displayed in Table 1. The optimum weight should decrease as the allowable control effort /3 2 is 
increased at constant allowable output response a 2 (left to right across the table), and the optimum weight 
should decrease as the allowable output response a 2 is increased at constant allowable control effort (3 2 (down 
the table). With reference to Table 1, we can see that this trend is observed in a macroscopic sense only, there 
being several examples where this trend is not observed. For example, considering the first column of Table 1, 
which corresponds to ft 1 = 50 lb 2 for varying a 2 , we see only two exceptions to the expected trends, these 
being at a 2 values of 6 x 10“ 5 and 9 x 10” 5 . Similar results are observed in all other columns and rows of 
Table 1. Results obtained using a state-space realization based on the physical displacements and velocities of 
nodal points are more consistent than when using the modal variables, although still not totally uniform. It 
might be pointed out that the results when using a physical realization were consistently easier to obtain, there 
being no need to alter the nominal value of j min to obtain convergence, and the number of global iterations 
required for convergence being consistenty lower. 

Some understanding of these contradictory results can be found by considering Table 2, which gives the 
optimal element areas found for /? 2 = 50 and for the varying a 2 corresponding to the first column of Table 1. 
Also given in this table is the number of global iterations required for convergence, the final values of the 
Lagrange multipliers (which then defines the LQR controller), and the initial value of the structural design 
variables (all the same for the symmetric set of initial conditions) at which the initial scaled system satisfies the 
constraints. Immediately apparent from Table 2 is a number of seemingly separate regions of the design space 
into which this structure has converged. For example, the final designs for a 2 = 5 x 10” 5 and a 2 = 7 x 10” 5 
seem to be similar in relative structure. Here, “similar” refers to the relative sizing of the structural members, 
in that design variables that are “larger” in one design are “larger” in the other. Both these designs are 
however distinctly different from those for a 2 = lx 10“ 5 and or 2 = 3 x 10~ 5 , which themselves are similar. 
The conclusion seems to be that we are converging into different regions of the design space with our solution 
algorithm, and that there are numerous local minima. Several columns of Table 2 seem to define their own 
region of the design space, being dissimilar to any other column. In other words, our design space seems to 
have multidimensional corrugations leading to multiple local minima. The solutions will lie somewhere on the 
intersection hyperplane between the surface of constant allowable output response and the surface of constant 
allowable control effort. 

This corrugated nature of the design space can be illustrated by considering the solutions obtained, for the 
same constraint case, when starting from different initial conditions. For the case of /? 2 = 75 and a 2 = 1 x 10~ 5 , 
Table 3 summarizes the results of runs made when modal state space realizations were used, and when only the 
initial conditions are varied. The different initial conditions are defined by setting all structural elements equal 
except the first (element 1), to which is added a percentage of the size of other elements. Even with this limited 
variation in the initial conditions, there are seemingly many distinct regions in the design space into which the 
system may converge. A picture of the constraint surfaces as a one-dimensional slice of the multidimensional 
space will emerge if these optimal structures are varied into each other in a linear fashion, and the constraint 
values are calculated between each case. That is, the structural design variables and Lagrange multipliers are 
changed linearly from the optimal values in one case to those in another case. Then the constraint surfaces 
obtained would be those seen when travelling in a straight line between each successive point. 

The results of such an analysis are shown in Figures 3 for the cases corresponding to those given in Table 3. 
As expected, the weight varies linearly between the cases, but it is the constraint curves that are much more 
revealing. For example, considering Figure 3, one can see that between case 1 and case 2, there is a “ridge” 
of output response larger than the maximum allowable value. Similarly, the control effort first decreases, then 
also increases to a ridge of high value. This corresponds to a hump in the constraint surfaces between the two 
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points in the design space. Assuming that we would see such behaviour when moving in every direction away 
from case 1 and case 2, rather than just in a direction between the two as shown in Figure 3, then the design 
points corresponding to these cases would represent local minima. In this situation, the design can become 
“trapped” in such a locally convex region, causing the solution algorithm to converge to different points. 

With reference to the same Figure 3, one can see that both the output response and control efforts are 
virtually constant between cases 2 and 3, while the weight increases slightly from 2053.0 lb to 2090.6 lb. 
This indicates that case 2 and case 3 actually represent the same optimal solution, with the difference being 
accounted for in the variance allowed by the convergence criteria used. The direction in the design space 
represented by the movement from case 2 to case 3 would lie in the intersection hyperplane of the surfaces of 
constant control effort and output response constraints, and would be at a shallow angle to the linear surface 
of constant weight. Figure 3 graphically illustrates a design space that is a very complicated function of the 
design variables, in which multiple local minima abound. 

There are some other tests that can be made on the hypothesis that the design is becoming trapped in 
local minima. If the design is actually trapped in a local minimum, the solution should stay in the vicinity of 
that minimum if the problem is changed only slightly. That is, if a converged solution is used as the initial 
conditions for an optimization run where the constraint objectives are changed by a “small” amount, then the 
new problem should converge to a point that is “close to” the initial point. Table 4 represents such a situation. 
Here, the solution was first obtained for the case where (3 2 = 50 and a 2 = 1 x 10” 5 , and where a modal state 
space realization and inverse design variable approximations were used. This converged solution was then used 
as the initial conditions for the cases f3 2 = 50 and a 2 = 2 x 10~ 5 , and f3 2 = 60 and a 2 = lx 10” 5 . Moving 
down each column, and across the top row, of Table 4, the converged solution from the previous case was used 
as the initial condition for the new problem. As can be seen from Table 4, the two expected trends in the data, 
as mentioned previously, are now observed without exception. The optimal solutions for the first column of 
Table 4, corresponding to the cases where /? 2 = 50, are given in Table 5 The solutions now appear to be in 
the same local region of the design space, as evidenced by the relative sizing of the optimal structures. For 
example, note that in all converged designs, structural elements 9, 10, and 12 are at their lower gage limit of 
0.1 in 2 , and that the first structural element is the largest by far. The optimal solutions for the /? 2 = 60 cases 
from Table 4 also appear to be in this same region of the design space. These results test the hypothesis that 
designs are converging to local minima, and indicate that the local optima are real and not simply figments of 
a numerical imagination. 

5,2 Direct Output Feedback Control 

Recall that when using direct output feedback, no simplifying assumptions can be made regarding the controller 
design variables (elements of K) ) such as the LQR assumption used in the case of full state feedback. Therefore, 
no scaling of the structure and controller to the closest constraint surface is performed. If a full state feedback 
case is to be solved in the form of direct output feedback with H — /, the number of design variables will 
increase significantly over the number of design variables created when other types of controllers are considered. 
This is because the number of states in the plant model will generally be large for anything but trivial systems, 
hence K will have many elements, all of which will be treated explicitly as design variables. However, such a 
situation is considered here to aid a comparison between the two solution algorithms, and the results obtained 
in the previous section. For the continuation algorithm, convergence was obtained after the specified number 
of global iterations, set by the specification of Ay. Also note that, since only one each of the output response 
and control effort constraints are specified at this stage (n a = 1, = 1), these can both be set as equality 

constraints without loss of generality since they will both be active at an optimum. 

Table 6 gives the optimal weights found using the sequential approximations solution algorithm in the case 
when a physical state space realization and inverse design variable approximations are used, for /3 2 = 50 and 
/ 3 2 = 60, and for varying a 2 . Compared to the similar case when scaling was performed, the optimal weights 
found here are larger in every case. Additionally, the solution times were significantly larger because of the 
number of iterations required for convergence. Note that some values in Table 6 are for situations where an 
average steady state value was obtained, but where the design was jumping around too much over each global 
iteration for convergence to occur. This was so even though the move-limits for every case in Table 6 were 
set at a relatively small ±2.5%. A smaller move-limit would aid convergence, but slow it considerably. Also, 
smaller move-limits may cause premature convergence if the design is at a point where the constraint surfaces 
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and the surface of constant weight are nearly parallel. 

Table 7 gives the optimal weights found using the continuation algorithm for the same cases listed in Table 6 . 
All cases listed were obtained with the continuation parameter Ay = 0.01 until y = 0.5, when A 7 became 0 . 02 , 
so that for every case convergence was achieved in 75 global iterations. Also listed in Table 7 are the move- 
limits on the maximum parameter changes allowable locally (e in Section 3.2). These are set to give a tradeoff 
between the satisfaction of the local linearity assumption and convergence to the local neigbouring problem, 
and must be found by numerical experimentation. Note that in both situations represented by Tables 6 and 
7, minimum move-limits on the elements of K are set so that these elements can change sign if desired. 

The solutions given in Tables 6 and 7 were obtained when the initial structure was defined with all truss 
elements of equal cross-sectional area. The particular values for this area are given in Tables 6 and 7 for 
each case, and were chosen so that the initial output response was "close to” its desired final value. With 
all structural elements at 120 in 2 , 90.0 in 2 , and 60.0 in 2 , the initial output responses were 2.067 x 10" 5 , 
3.674 x 10- 5 , and 8.266 x 10 “ 5 respectively. The initial control effort was zero of course, since all elements of 
K were set initially to zero. 

Comparing the optimal weights from Tables 6 and 7, it can be seen that those obtained using the con- 
tinuation solution algorithm are significantly lower than those obtained using the sequential approximations 
solution algorithm in every case. The optimal weights are also much more consistent, in terms of the expected 
trends as a 2 increases, when using the continuation method. Convergence was obtained for every case listed 
in Table 7, whereas for three of the cases listed in Table 6 , no convergence was obtained within 300 iterations. 
The reason for the convergence failure can be illustrated by considering the convergence histories given in 
Figures 4 and 5, for typical cases from Tables 6 and 7 respectively. The inset in Figure 4 shows some of the 
histories toward the end of the solution in more detail. These can be seen to be very rough, as compared to 
the histories in Figure 5 which are smooth everywhere. These parameter oscillations for the sequential approx- 
imations solution technique are indicative of a move-limit set too high, so that the local linearity assumption 
is violated. However, since the convergence is very flat toward the end of the solutions, a smaller move-limit is 
likely to cause premature convergence, or to significantly slow down the convergence for very little additional 
objective reduction. 

The optimal weights found using the continuation method listed in Table 7 compare very favourably with 
those found for the same cases (/? 2 = 50 and varying a 2 ) when the full state feedback LQR assumption was 
used to simplify solution. In almost every case, the optimal weight is approximately equal to or lower than 
those found earlier. However, the sequential approximations solution algorithm results are much worse, being 
significantly larger than the optimal weights found earlier in every case. The sequential approximations solution 
algorithm performs so poorly because it tends to become trapped in a local minimum close to the specified 
initial conditions. 

Consider now cases where the full state is not available for feedback. In the DRAPER I structure, the six 
right-angled bipods are usually assumed to take on the duties of force actuators and colocated rate sensors, so 
that H = G 7 (output state of dimension six). Tables 8 and 9 list the optimal weights found for this reduced- 
order output vector for the same cases used in Table 6 , when the sequential approximations and continuation 
solution algorithms respectively are used. Note that the optimal weights consistently obey the trends that are 
expected as the allowable output response and controller effort are altered. Even so, the designs can converge 
into completely different regions of the design space for any particular case. This is illustrated by Table 10, 
which gives the optimal structural design variables for the cases represented by the second column of Table 9. 
Note that every solution does not define a unique local minima. Many of the solutions seem to lie in the same 
local minima region, for example the solution for the case a 2 = 1,7, and 9x 10~ 5 . Once again, the design space 
is found to have many local minima. 

The two solution algorithms here seem to predict quite similar optimal weights for the cases considered, 
although in general those found using the continuation method are slightly better. However, convergence for 
the results using the sequential approximations solution algorithm were more difficult to obtain than those for 
the continuation method. Altering the nominal move limits (set at ±2.5%), and perhaps reducing it toward 
the latter stages of each solution would aid convergence, but at increased effort on the part of the user. The 
continuation results are easier to obtain since they require less individualized attention. 

Comparing Tables 8 and 9 to Table 7, it can be seen that the optimal weights found when H — G T are 
much larger than those found when full state feedback was used. The reason for this is the particular placement 
of the disturbance forces relative to the design output states (which determine the output response). For the 
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case considered here, the DRAPER I structure is disturbed at node 1, the displacement of which is to be 
kept below the design objective value a 2 . With full state feedback, the displacement states of node 1 are 
available for feedback, whereas if H = G T , these states are not available, and the effect of the displacement 
of node 1 is available only indirectly through its effect on the velocities along the six bipods that make up 
the sensors. The importance of these states can be seen by examining the optimal gain matrices for the full 
state feedback case. The largest gains associated with displacements and velocities in these matrices appear 
in the columns corresponding to node 1 degrees of freedom, indicating that the states associated with node 1 
are very important. When they are not available, the controller does not have as much information about the 
state of node 1, the node it is trying to control, as it does in the case of full state feedback, and will increase 
the structural stiffness (and hence mass) to compensate. 

To illustrate further, consider output feedback with the displacement and velocity states of node 1 added 
to the output vector used previously. The optimal weights for the same (3 2 = 50 and f3 2 = 60 cases used in 
Table 9 with this new output vector (of dimension 12 now) are given in Tables 11 and 12, when the sequential 
approximations and continuation solution algorithms are used respectively. There are now larger differences 
between the weights found using the two solution algorithms than in the case when H = G? only, with the 
continuation method giving the best results (with one exception). The weights obtained using the continuation 
algorithm are now very close to those obtained when using full state feedback (in Table 7), although still a 
little larger in every case. This is because even with this larger output vector, the displacement states for 
nodes two through four are still not used in the controller. 

All the results presented so far were generated with the external stochastic disturbance intensities set at 
one ( X w = 7). If these intensities are varied as X w = x w I , the effect of varying x w on the optimal weight is 
shown in Figure 6. These results were generated using the continuation solution algorithm, for the case when 
a 2 = 1 x 10“ 5 and /3 2 = 50, and when all structural design variables were initially set at 120 in 2 . As can be 
seen, the relationship between the disturbance intensity and the optimal structural weight appears basically 
linear in the range shown. However, a linear fit of the data does not produce an optimum weight of zero for 
a zero intensity disturbance, as would be expected. Therefore, the relationship cannot be exactly linear. The 
results also indicate that the disturbance level chosen for the previous results (a^ = 1) was significant for the 
range of output response and control effort constraint objectives used. 

As long as the neighbouring problems in the continuation algorithm were close enough, which was satisfied 
by starting from an initial point where the constraint values were “close to” their final desired locations, no 
difficulties were experienced in obtaining convergence. The solution times for the continuation algorithm were 
2-3 times longer than for solutions by the sequential approximations algorithm, since on average approximately 
8-12 local iterations were required for convergence to the neighbouring problem for each global iteration. The 
performance of the sequential approximations solution algorithm decreases as the dimensionality increases, as 
evidenced by the sequence of cases where H — G? (48 design variables), H = G T plus node 1 displacement 
and velocity states (84 design variables), and H = I (156 design variables). However, the continuation method 
seemed much less sensitive to the problem dimension. The continuation solution method was found to be 
generally superior, for this problem at least, to the sequential approximations solution method, with respect to 
the confidence in obtaining a “good” converged solution. The disparity in solution times was acceptable because 
of the ease with which solutions were obtained using the continuation method, and because the solutions found 
seemed to be generally much better than those found using the sequential approximations algorithm. 


6 Conclusions 

In this work, the integrated control /structure design optimization problem has been investigated from a re- 
sponse to disturbances point of view. Both full state and output feedback controllers were employed in the 
control strategy, and two solution methods were compared. It was found that for this problem, the continuation 
method coupled with a sequential linear programming approach performed better than the more traditional 
type of nonlinear approximations approach, in the sense that it was more robust to changes in the arbitrary 
parameters set by the user, obtained better results, and the results were easier to obtain. The design space was 
found to exhibit multiple local minima in which the solution could become trapped, although the continuation 
solution method seemed to handle the corrugated design space better than the other method. In future work, 
more diverse controller types must be considered, along with structures consisting of more complicated finite 


134 



elements than simple truss members. Additionally, more realistic problems of higher dimension must be solved, 
to demonstrate the practicality of this design procedure. 
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Figure 1: The DRAPER I Structure 



Figure 2: The disturbance model 
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5 
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6 

1460.6 
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1009.8 

7 

1209.0 

1076.0 

928.7 

1052.7 

8 

1066.7 

1214.1 

987.9 
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9 

1126.3 

940.5 
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720.1 

10 

971.9 

1035.1 

766.4 

970.6 


Table 1: Optimal weight using a modal stale-space realization, the 
symmetric set of initial conditions, and inverse design variable approx- 
imations. 
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